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Abstract: One of the unknown parameters in neutrino oscillation studies is the octant of 
the atmospheric neutrino mixing angle 023- In this paper, we discuss the possibility of de- 
r"! ■ termining the octant of 823 in the long baseline experiments T2K and NOz^A in conjunction 

with future atmospheric neutrino detectors, in the light of non-zero value of ^13 measured 
by reactor experiments. We consider two detector technologies for atmospheric neutrinos - 
magnetized iron calorimeter and non-magnetized Liquid Argon Time Projection Chamber. 



> 
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I We present the octant sensitivity for T2K/N0z^A and atmospheric neutrino experiments 

separately as well as the combined sensitivity. For the long baseline experiments, a precise 
' measurement of ^13, which can exclude degenerate solutions in the wrong octant, increases 

■ the sensitivity drastically. For 623 = 39° and sin^ 2^13 = 0.1, at least ~ 2a sensitivity can 

be achieved by T2K + NOj^A for all values of 5cp for both normal and inverted hierar- 
chy. For atmospheric neutrinos, the moderately large value of ^13 measured in the reactor 
experiments is conducive to octant sensitivity because of enhanced matter effects. A mag- 
netized iron detector can give a 2a octant sensitivity for 500 kT yr exposure for 623 = 39°, 
Scp = and normal hierarchy. This increases to 3a for both hierarchies by combining 
with T2K and NOz^A. This is due to a preference of different 623 values at the minimum 
by T2K/N0z^A and atmospheric neutrino experiments. A Liquid Argon type detector 
for atmospheric neutrinos with the same exposure can give higher octant sensitivity, due 
to the interplay of muon and electron contributions and superior resolutions. We obtain a 
~ 3a sensitivity for 623 = 39° for normal hierarchy. This increases to ^ 4(T for all values 
of 6c p if combined with T2K/N0z^A. For inverted hierarchy the combined sensitivity is 
around 3a. 
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1. Introduction 

The measurement of a non-zero value of the mixing angle ^13 by the reactor experiments 
Double-Chooz [1], Daya-Bay [2] and RENO [3] heralds a major breakthrough in the ad- 
vancement of neutrino physics. The best-fit value of sin^ ^13 from latest global analysis of 
solar, atmospheric, reactor and accelerator data is 0.023 ± 0.0023 [4]. which signifies non- 
zero ^13 at 10(7 level. Other global analyses also give similar results [5,6]. This confirms 
the earlier observation of non-zero ^13 in T2K [7] and MINOS [8] experiments as well the 
indication of non-zero best-fit value of ^13 from previous analyses [9,10]. Among the other 
oscillation parameters the solar parameters Am2i and 612 are already well measured from 
solar neutrino and KamLAND experiments [11,12]. The solar matter effect also dictates 
Arrigi > 0. The most stringent constraint on mass squared difference Am|^ governing the 
atmospheric neutrino oscillations comes from the data from the MINOS experiment [13]. 
However the ordering of the third mass eigenstate with respect to the other two i.e. the 
sign of Amg]^ is not yet known. There are two possible arrangements of the neutrino mass 
states: (i) mi < m2 < corresponding to Normal Hierarchy (NH) and (ii) < mi < m2 
corresponding to Inverted Hierarchy (IH). The mixing angle sin^ 622, is mainly determined 
by the SuperKamiokande (SK) atmospheric neutrino data. However the octant in which 
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parameter 


best fit 


3cr range 


Amii [10-5 eV^] 


7.54 


6.99-8.18 


[Am^J [10-3 eV^] 


2.43 
2.42 


2.19 - 2.62 
2.17-2.61 


sin^ 012 


0.307 


0.26-0.36 




0.386 
0.392 


0.33-0.64 
0.34-0.66 


sin^ 6^13 


0.0241 
0.0244 


0.017-0.031 


S 


l.OSvr 
l.ODvr 


- 27r 



Table 1: The best-fit values and Str ranges of neutrino oscillation parameters from Ref. [5]. 
The first and second rows in each case denote the values for a normal and inverted neutrino mass 
hierarchy respectively. 

this mixing angle lies is not yet decisively determined by the data. A full three-flavour fit to 
the SK data gives the best-fit for NH in the lower octant (LO) and IH in the higher octant 
(HO) keeping ^13 as a free parameter in the analysis [14]. Three-flavour global analysis 
of all available neutrino data give the best-fit 623 in the lower octant. The best-fit values 
and 3(7 ranges of oscillation parameters from global analysis in Ref. [5] is summarized in 
Table |l|. At present there is no significant constraint on the CP phase 6c p and the whole 
range from — 27r is allowed at the 2a level. 

There are two aspects in the measurement of ^13 which is expected to play an important 
role in resolving the outstanding issues in neutrino oscillation physics. Firstly if ^13 is non- 
zero and relatively large — (i) CP violation in the lepton sector can be probed (ii) the 
earth matter effect on the propagation of neutrinos can be sizable. The latter facilitates 
the determination of mass hierarchy and octant in experiments in which neutrinos travel 
through an appreciable path length. The second aspect is the precision measurement of 
013 : which helps in increasing the sensitivity in the determination of hierarchy, octant and 
6c P- 

The octant degeneracy means impossibility of distinguishing between 023 and 7r/2 — 023- 
This is generic and robust for vacuum oscillation probabilities that are functions of sin'^ 2023 
e.g. the two flavour muon survival probability in vacuum [15]. If on the other hand the 
leading term in the probability are functions of sin^ 023 (e.g. -P/^e) then the inherent octant 
degeneracy is not there but lack of knowledge of other parameters like 0i3 and 6c p can give 
rise to octant degeneracy [16, 17]. These issues may affect the octant sensitivity of the long 
baseline experiments T2K and NOz^A where the matter effect is not very significant and in 
particular resonant matter effects do not get a chance to develop. Although conventionally 
the octant degeneracy refers to the indistinguishability of 023 and 7r/2 — 023, in view of the 
present uncertainty in the measurement of 023 the scope of this can be generalized to any 
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value of 023 in the wrong octant within its allowed range. In this paper by octant sensitivity 
we refer to this "generalized" definition. If in addition the hierarchy is unknown then there 
can also be the wrong-octant-wrong-hierarchy solutions and one needs to marginalize over 
the wrong hierarchy as well. 

Atmospheric neutrinos pass through long distances in matter and they span a wide 
range in energy and can encounter resonant matter effects. In this case the octant sensitivity 
in P^^ ensues from the term sin^ 613 sin^ 20^ [18]. Pfie in matter contains sin^ ^13 sin^ 2^^. 
Since at resonance sin^ 2^^^ « 1, the octant degeneracy can be removed. In this case also 
one can probe the effect of 6cp uncertainty on the lifting of this degeneracy. 

Atmospheric neutrinos provide fluxes of both neutrinos and antineutrinos as well as 
neutrinos of both electron and muon flavour. On one hand it provides the advantage of 
observing both electron and muon events. However on the other hand a particular type of 
event gets contributions from both disappearance and appearance probabilities. This can 
be a problem if the matter effects for these two channels go in opposite directions. Thus 
it is necessary to carefully study the various contributions and ascertain what may be the 
best possibility to decode the imprint of matter effects in atmospheric neutrino propaga- 
tion. Three major types of detector technologies are under consideration at the present 
moment as future detector of atmospheric neutrinos. 

(i) Water Ccrcnkov detectors : Such type of detectors have already been shown to be a 
successful option for atmospheric neutrino detection by the SuperKamiokande collabora- 
tion. This is sensitive to both electron and muon events and the energy threshold can 
be relatively low. Megaton detectors of this kind under consideration for future are HK, 
MEMPHYS [19,20]. These cannot be magnetized and hence provide no charge identifica- 
tion capability. Multi-megaton detectors with ice also fall in this category. An example of 
such a detector is PINGU [21,22], which is a proposed upgrade of the DeepCore section of 
the IceCube detector for ultra-high energy neutrinos [23] with a lower energy threshold for 
atmospheric neutrino detection. 

(ii) Magnetized Iron Detectors: Such a detector for atmospheric neutrinos were proposed 
by the MONOLITH [24] collaboration and is now actively pursued by the INO collabora- 
tion [25]. This has a relatively high threshold and is mainly sensitive to muon neutrinos. 
These type of detectors offer the possibility of magnetization, thus making it possible to 
distinguish between muon and antimuon events, which enhances the sensitivity. 

(iii) Liquid Argon Time Projection Chamber (LArTPC): Examples of such detectors are 
ICARUS and ArgoNeuT [26,27]. The hallmark of these detectors are their superior particle 
identification capability and excellent energy and angular resolution. They are sensitive to 
both electron and muon events with good energy and direction reconstruction capacity for 
both type of events [28]. Possibility of magnetization is also being discussed [29,30]. 

In this paper we study in detail the possibility of removal of octant degeneracy in 
view of the precise measurement of a relatively large value of ^13. There have been many 
earlier studies dealing with this subject both in the context of long baseline and atmo- 
spheric neutrinos. The importance of combining accelerator and reactor experiments and 
the role of a precise measurement of 613 in abating the octant degeneracy have been con- 
sidered in Refs. [31-33]. The possibility of resolving the octant degeneracy using long 
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baseline experiments has also been explored in Refs. [34-39]. Recently octant sensitivity 
in the T2K/N0i/A experiments has been investigated including the recent results on the 
measurement of non-zero ^13 by reactor experiments [40] . The octant sensitivity for atmo- 
spheric neutrinos in the context of magnetized iron calorimeter detectors was considered 
in Refs. [18,41,42], for water Cerenkov detectors in Ref. [19] and for LArTPC in Ref. [43]. 

We examine the octant sensitivity in the long baseline experiments T2K and NOz^A 
and in the atmospheric neutrino experiments as well as the combined sensitivity of these 
experiments. In particular we address whether degeneracy due to ^13 can still affect octant 
determination at the current level of precision of this parameter. For fixed values of ^13 
the effect of lack of knowledge of 5c p on the octant determination capability of these 
experiments is also studied. We take into account the uncertainty of 622, in the wrong 
octant and discuss how much this can influence the octant sensitivity. We present results 
for the two cases of known and unknown hierarchy. 

For the study of atmospheric neutrinos we consider magnetized iron calorimeter detec- 
tors with charge sensitivity which is sensitive to the muon neutrinos. We also consider a 
non-magnetized LArTPC detector which can detect both electron and muon neutrinos. In 
particular we discuss the interplay between the muon and electron type events in the overall 
octant sensitivity. For our atmospheric analysis we assume a prior knowledge of hierarchy. 
Lastly we do a combined analysis of T2K, NOz^A and atmospheric neutrinos and discuss 
the synergistic aspects between long baseline and atmospheric neutrino experiments. 

The plan of the paper is as follows. In Section ^ we discuss the octant degeneracy at 
the level of oscillation and survival probabilities, for baselines corresponding to both atmo- 
spheric neutrinos and those relevant to NOz^A and T2K. Section ^ discusses the analysis 
procedure and results. First we discuss the octant sensitivity in NOj^A and T2K. Next we 
describe the results obtained for octant sensitivity using atmospheric neutrino detectors. 
Finally we present combined octant sensitivity of long baseline and atmospheric neutrino 
experiments. We end by summarizing the results. 

2. Analysis of octant degeneracy 

The ambiguity in the determination of the octant of ^23 niay appear in the oscillation and 
survival probabilities as 

(a) the intrinsic octant degeneracy, in which the probability is a function of sin^ 2023 
and hence the measurement cannot distinguish between ^23 and 7r/2 — ^23; 

p{e%) = p{^/2-e%) (2.1) 

(b) the degeneracy of the octant with other neutrino parameters, which confuses octant 
determination due to the uncertainty in these parameters. In particular, this degeneracy 
arises in probabilities that are functions of sin'^ ^23 or cos^023- For such cases P{02s) 7^ 
P{tt/2 — ^Is)- However for different values of the parameters ^13 and 5c p the probability 
functions become identical for values of ^23 in opposite octants for different values of these 
parameters, i.e. 

P(.0il,ei3,5cp) = P(vr/2 - e%,9[„5'cp), (2.2) 
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where denotes the true value of the mixing angle and the primed and unprimed values 
of and 5c p lie within the current allowed ranges of these parameters. In the case of 
5cp, this covers the entire range from to 2it, while for 613 the current 3a range is given 
by sin^ 2^13 = 0.07 — 0.13 . From the above equation it is evident that even if ^13 is 
determined very precisely, this degeneracy can still remain due to complete uncertainty in 
the CP phase. In fact, the scope of this degeneracy can be enlarged to define this as 

P{eil, e,s, Sep) = P{e^r'^ ^'13, ^cp) (2.3) 

where 62^°^^ denote values of the mixing angle in the opposite octant. 

The features of the octant degeneracy and the potential for its resolution in differ- 
ent neutrino energy and baseline ranges can be understood from the expressions for the 
oscillation and survival probabilities relevant to specific ranges. We discuss below the 
probability expressions in the context of the fixed baseline experiments N0i/A/T2K and 
for atmospheric neutrino experiments. 



2.1 Neutrino Propagation in Matter 

Neutrinos travelling through earth encounter a potential due to matter given as, 



A = 2\/2GFncE = 2x 0.76 x lO"'' x Yc 



P 



y/cc 



E 



GeV 



eV^ 



(2.4) 



where Gp is the Fermi coupling constant and rig is the electron number density in matter, 
given by n^; = NaY^p {Na = Avogadro's number, Yg = electron fraction ~ 0.5, p = earth 
matter density. 

The mass squared difference (Am|^)™ and mixing angle sin^ 2^5^ in matter are related 
to their vacuum values by 



(Am2^ cos 26li3 - A)2 + (Am^^ sin 26*13)2 



sin 26* 



13 



Ato^i sin 2^13 



^(Am2^ cos 2^13 - A)2 + (Am^^ sin 2^13)2 



(2.5) 



The MSW matter resonance [44-46] occurs and the mixing angle 0"^ becomes max- 
imal at neutrino energies and baselines for which the terms Amg^^ cos 2^13 and A in the 



denominator of eq.(2^) become equal. Hence the matter resonance energy Ej-es is given by 



|Am2^| cos 29 



13 



2 X 0.76 X 10-4 X Yc X p 



(2.6) 



Since the corresponding expression for antineutrinos is obtained by making the re- 
placement A — )• —A, it may be observed that the matter resonance occurs for the normal 
neutrino mass hierarchy (i.e. Am|^ > 0) for neutrinos and for the inverted mass hierarchy 
(i.e. Am|^ < 0) for antineutrinos. 
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2.2 Octant ambiguity in P^g and P^^ 



For N0z/A/T2K, the baselines are 812 and 295 Km respectively and the peak energies of the 
beams are in the range 0.5-2 GeV. For these values of baselines, the earth matter density 
is in the range 2.3 — 2.5 g/cc, and the corresponding matter resonance energies are above 

10 GeV. Hence the neutrino energies of both experiments lie well below matter resonance, 
and the oscillation probabilities will only display small sub-leading matter effects. The 
expressions of the relevant probabilities and P^e in vacuum are given by the following 
expressions obtained in the one-mass scale dominant (OMSD) approximation. 



pv 



1 - sin^ 26*23 sin^ 



1.27 Amii I 



+ 4sin^6'i3 sin^6'23 cos 2^23 sin^ 



1.27 Am^, I 



(2.7) 



sin ^23 sin 2^13 sin 



1.27 Am^i - 



E 



(2.8) 



The above probability expressions have sub-leading corrections corresponding to small 
matter effect terms and the solar mass-squared difference [47-51]. We observe the following 
salient features from these expressions: 

(a) The disappearance channel has in its leading order a dependence on sin^ 2^23, 
and hence is dominated by the intrinsic octant degeneracy. There is a small ^13-dependent 
correction in the measurement of ^23 which gives a minor (~ 1%) resolution of the degen- 
eracy if 9i3 is known precisely. 

(b) The appearance channel P^e has the combination of parameters sin^ ^23 sin^ 2^13, 

and hence docs not STiffcr from the intrinsic octant degeneracy. However, the degeneracy 
of the octant with the parameter ^13 comes into play, since the above combination may be 
invariant for opposite octants for different values of ^13, and hence this degeneracy cannot 
get lifted with a measurement from such experiments alone [32]. This channel can be also 
affected by the large uncertainty in Scp when sub-leading corrections are included. 

For atmospheric neutrinos, the relevant baselines and energies are in the range 1000 
- 12500 Km and 1-10 GeV respectively. A large region in this L and E space exhibits 
strong resonant earth matter effects, since the earth densities in this baseline range (3-8 
g/cc) correspond to resonance energies Efcs = 4 — 9 GeV. Hence the relevant probability 
expressions P™ , P™ and PJ^ can be written, in the OMSD approximation and with full 
matter effects, as [48,52] 



P;7e = sin^ 023 sin^ 26^^ sin^ 



1.27 (Amii) 



E 



(2.9) 
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- cos2 0- sin2 2023 sin^ ^1.27 

- sin2 0?^3sm2 2023 sin^ 1.27 (^^I^l±l 

- sin^023 sin^ 205^3 sin^ 1.27 (Am^i)"^ ^ 




E 



E 



L 



E 



L 



(2.10) 



.m 

ee 



2 205^3 sin^ 1-27 (Am^^) 



m 



L 



(2.11) 



1 — sin' 



E 



In this case, the following features are observed: 

(a) The oscillation probability in matter P™ is still guided to leading order by a 
dependence on sin^023- But strong resonant earth matter effects help in resolving the 
degeneracy since the mixing angle ^13 in matter gets amplified to maximal values (close 
to 45'') near resonance. The combination sin^ ^23 sin^ 20^^ no longer remains invariant over 
opposite octants, since sin^ 20^ becomes close to 1 in both octants irrespective of the 
vacuum value of ^13. This breaks the degeneracy of the octant with ^13. 

(b) The muon survival probability in matter P™^ has leading terms proportional to 
sin^2023) as in the vacuum case, which could give rise to the intrinsic octant degeneracy. 
But the strong octant-sensitive behaviour of the term sin^ ^23 sin^ '^^13 near resonance can 
override the degeneracy present in the sin^ 2023-dependent terms. 

(c) The electron survival probability P™ is independent of ^23 and hence does not 
contribute to the octant sensitivity. 

(d) Since the P™ channel is simply the CP conjugate of P™ , the probability level 
discussion in this section is applicable to P™ also. 

In the following discussion, we address the octant degeneracy due to 023 in the wrong 
octant, and unknown values of dcp, at a probability level. The probability figures |l| - 
are drawn by solving the full three flavour propagation equation of the neutrinos in matter 
using PREM density profile [53]. In all these figures the left panels are for the NOi^A peak 
energy and baseline (2 GeV, 812 Km), while the right panels are for a typical atmospheric 
neutrino energy and baseline (6 GeV, 5000 Km). The top row denotes the appearance 
channel P^e, while the bottom row denotes the disappearance channel i^^/^- 

Figure |l] depicts the probabilities P^e and P^^ as a function of sin^ 26*13 for 9^^^ = 39° 
and 62^"""^ = 51°. The bands show the probability range in each case when 6c p is varied 
over its full range (0 to 27r). For a given fixed value of sin'^ 20i3, the distinction between 
023 in the two octants can be gauged from the separation of the two bands along the 
relevant vertical line. The left panels show that for P^e, the 6c p bands overlap and the 
two hierarchies cannot be distinguished till nearly sin^ 20i3 = 0.1. The inset in the upper 
left panel shows the region of separation of the bands near sin^ 20i3 = 0.1 in detail. Hence 
the knowledge of the parameter 0i3 upto its current level of precision becomes crucial, 
since a 0i3 range including lower values would wash out the octant sensitivity derivable 
from such experiments due to the combined degeneracy with 0i3 and 6c p- For P^^, the 
intrinsic degeneracy predominates and the effect of 6cp variation is insignificant. 
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Figure 1: Behaviour of tlie muon survival and oscillation probabilities as a function of sin^ 20i3 
showing the 023 octant degeneracy and its breaking. The left panels are for the NOi/A peak energy 
and baseline (2 GeV, 812 Km), while the right panels are for a typical atmospheric neutrino energy 
and baseline (6 GeV, 5000 Km). The top row denotes the appearance channels P^e and Pg^, while 
the bottom row denotes the disappearance channel . The values of oscillation parameters chosen 
are = 39°, 623°^^ = 51°. The bands denote a variation over the full range of the phase Sep- 
The inset shows the region of separation of the bands near sin^ 2^13 =0.1. 



For the 5000 km baseline, due to the resonant matter effects breaking the octant 
degeneracy at the leading order, both P^e and P^^ show a wider separation between the 
opposite-octant bands, even for small values of 613. This is due to the sin^ ^23 sin^ 20^^ 
(sin'* 023 sin^ 20]^) term in P^e (P^t/.t)- The Scp bands in the right-hand panels are much 
wider because of the enhancement of the subleading terms due to matter effects. However, 
the enhancement is more for the leading order term which alleviates the degeneracy with 
Scp- 

Figure |2| shows the probabilities P^e and P^^ as a function of sin^ 2^13 . We have held 
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Figure 2: Same as Figure ^withi a fixed true and test Sep — and true 6*13 = 39°, witli tire band 
denoting a variation over tlie full allowed range of 6*^3'^°"^ = 45° to 54° in the wrong octant. 



fixed, true and test 5cp = and true ^23 = 39°, with the band denoting a variation 
over the full allowed range of 62^°^^ = 45° to 54° in the wrong octant. Thus this figure 
reveals the effect of the uncertainty in the measurement of ^23 in the determination of 
octant, for a given value of sin^ 2^13. For the NOz^A baseline, the survival probability P^^ 
shows an overlap of the test 622, band with the true curve, while in the probability P^e 
there is a small separation. Figures |l] and |2| thus indicate that for the NOi^A baseline, the 
octant sensitivity from P^e is more affected by the uncertainty in 5c p and less by the test 
023 variation, while for P^^ the opposite is true. The plots for the atmospheric neutrino 
baseline show a clear breaking of the octant degeneracy in both P^o and P^^ even for small 
values of ^13 , indicating that the octant sensitivity from the atmospheric neutrino signal is 
stable against the variation of both 5c p and the test value of ^23; even for small values of 

^13- 
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Figure 3: Energy spectra of the probabilities P^c and P^^ for the NOi^A basehne (left panels) and 
for a sample atmospheric neutrino baseline 5000 Km (right panels) . The figure shows the variation 
in the probabilities in each case when sin^ 2^13 is varied over three values in the current allowed 
range, fixing the true and test 6*23 values 6*23 — 39° and 623°"^ = 51°. Sep is fixed to 0. 

Having discussed octant sensitivity in the context of fixed ^13, we now come to the 
question of uncertainty in this parameter. The current reactor measurements of sin^ 2^13 
have measured this parameter with a precision of 0.01 and hence the effect of ^13 uncertainty 
on octant degeneracy is largely reduced. In Figure |, the energy spectra of the probabilities 
P^e and P^^ are plotted. The figure shows the variation in the probabilities in each case 
when sin^ 2^13 is varied over three values - sin^ 2^13 = 0.07,0.1 and 0.13 covering the 
current allowed range, and for two illustrative values of O23 - 39° and 51° in opposite 
octants. The value of 6cp is fixed to be 0. The figure indicates that the separation of the 
LO and HO curves and hence the octant sensitivity depends on the value of sin^ 2^13 in 
opposite ways depending on whether the true 623 value lies in the higher or lower octant. 
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E (GeV) E (GeV) 

Figure 4: Energy spectra of the probabilities P^c and P^^ for the NOj^A basehne (left panels) and 
for a sample atmospheric neutrino baseline 5000 Km (right panels) . The figure shows the variation 
in the probabilities in each case when 9^^°"'^ is varied over the entire allowed range in the wrong 
octant, as well as varying the test Sep, fixing the true values = 39° and S^p = in the solid 
curve. The spread due to variation of both test parameters is denoted by the blue band. The black 
band shows the variation with test Sep for a fixed 02^°"^ = 51°. 

In P^e, for = 39°, lower values of true sin^ ^13 can give a higher octant sensitivity 
since they are more separated from the band of variation of the probability over the whole 
range of test sin'^ 2^13 , and as true sin^ 2^13 increases the degeneracy with the wrong 623 
band becomes more prominent. For ^23 = 51°, the opposite is true, i.e. higher values of 
sin^ 20^3 have a better separation with the wrong 623 band. For example, for the NOi^A 
baseline, if = 39°, the 

P^e curve for sin 26^^ — 0.07 is well-separated from the band 
for 0^^°"^ = 51, while the curve for sin^ 20^3 = 0.1 is seen to be just separated from it, 
and the sin^ 2,6^13 ~ ^-^^ curve lies entirely within the band of test sin^ 2^13 variation with 
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the wrong octant. But if 6*23 = 51°, the sin^ 20^3 = 0.07 curve suffers from the degeneracy 
while the curves for sin^ 20I3 = 0.1 and upwards he clearly outside the 62^"^^ band. For 
P^^, the effect of ^13 on the separation between the opposite octant bands is less, since the 
behaviour is governed by the intrinsic octant degeneracy. 

For the 5000 km baseline, due to strong matter effects, the separation of the true and 
wrong ^23 bands is much better, and only the highest (lowest) values of sin^ 20^3 suffer 
from a degeneracy with the wrong octant for P^e if the true octant is lower (higher). For 
P^^, the behaviour is reversed. 

Figure § again depicts the energy spectra of the probabilities P^e and P^^. This figure 
shows the variation in the probabilities in each case when 62^°^^ is varied over the entire 
allowed range in the wrong octant, as well as varying the test 5c p-, fixing = 39° and 
5qp = 0. The solid black band denotes the variation with test 5cp for a fixed 02^"^^ = 51°. 

For the NOi^A baseline, the test probability bands show an almost complete overlap 
with the ^23 = 39°,(5^p = curve. However, the octant sensitivity may still be retained 
due to spectral information. So there always exist specific energy regions and bins from 
which the octant sensitivity can be derived. 

At 5000 km, the minimum separation between octants does not occur at or near 
0wrong _ _ ^tr _ g-|^o seen from the solid black shaded region. The edge of 

the striped blue band corresponding to some other value of 623 is closest to the 'true' curve 
and well-separated from it. This shows that there is octant sensitivity for this baseline 
even after including the uncertainty in 5c p and 023- The 62^°"^ for which the minimum 
separation is likely to occur will be further discussed in the context of Figure |9|. 



3. Analysis and Results 

In this section we present the results of our analysis for T2K, NOz^A and atmospheric 
experiments. We also give results for octant sensitivity when the results from both type of 
experiments are combined. 

In Refs. [54,55], it has been shown that the atmospheric parameters Am^^^ and 
measured in MINOS are related to the oscillation parameters in nature, Am|^ and 623 
using the following non-trivial transformations: 

sm6'23 — TT- , (o-lj 

COS di3 

Amg]^ = Am^^^ + (cos^ 612 — cos 5 sin ^13 sin 2^12 tan ^23) Amg^ . (3.2) 

These transformations become significant in light of the moderately large measured value of 
613. Therefore, in order to avoid getting an erroneous estimate of octant sensitivity, we take 
these 'corrected' definitions into account. Thus, in calculating oscillation probabilities, we 
use the corrected parameters A.m'^i and 623, after allocating the measured values to Am^^^ 
and 9fj,^. 

Our analysis procedure consists of simulating the experimental data for some specific 
values of the known oscillation parameters known as the 'true values'. The experimental 
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data is generated for a fixed hierarchy and for the foUowing fixed values of the parameters: 



(Amii)*'^ = 7.6 X 10-5 eV^ 
(sin^ 012)*'^ = 0.31 
(sin^ 2^13)*'^ = 0.1 (3.3) 
{AmltJ'^ = 2.4 X 10-=^ eV^ , 

and specific values of Oj^^^ and S^p. In the theoretical predictions which are fitted to the 
simulated experimental data, the 'test' parameters are marginalized over the following 
ranges: 

6cp G [0,27r) 

(35°, 45°) (true higher octant) 
(45°, 55°) (true lower octant) 

sin^26'i3 G (0.07,0.13) . 

(3.4) 

Am|^ and sin^ 612 are fixed to their true values since their effect is negligible. Also, after 
verifying that the effect of a marginalization over Am^^ is minimal, we have fixed Am^^^ 
to its true value for computational convenience. 

In our calculation, the hierarchy is assumed to be known in all the atmospheric neutrino 
experiments, since the time-scale involved would ensure that the hierarchy is determined 
before any significant octant sensitivity is achievable. For N0i^A/T2K, results are given 
both with and without prior knowledge of hierarchy, i.e. marginalizing over the test hier- 
archy. 

Priors are taken in terms of the measured quantities sin^ 2^13 and sin^ 29^^ as follows: 



2 _ . sin^2g^7-sin^2g^^ \ / sin^ 2^1^ " sin^ 2^13 . 



with the la error ranges as crsin-^20i^^ ~ '^sm2 26»i3 ~ 0.01. The latter is the error on 

^13 quoted recently by Double Chooz, Daya Bay and RENO [1-3]. 

3.1 NOi/A and T2K 

We simulate the current generation of long baseline experiments NOz/A and T2K, using the 
GLoBES package [56,57] and its associated data files [58,59]. For T2K, we assume a 3 year 
run with neutrinos alone, running with beam power of 0.77 MW throughout. (We choose 
this low running time for T2K to compensate for the fact that their beam power will be 
increased to its proposed value over a period of a few years [60].) The energy resolutions 
and backgrounds are taken from Refs. [61-65]. The detector mass is taken to be 22.5 kT. 
For NOz^A, we consider the set-up as described in Refs. [66,67], which is re-optimized for 
the moderately large measured value of ^13. The 14 kT detector receives a neutrino and 
antineutrino beam for 3 years each, from the NuMI beam. 
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Figure 5: Marginalized octant sensitivity from a combination of NOt^A and T2K, for the case of 
normal and inverted mass hierarchy. The test hierarchy has been fixed to be the same as the true 
hierarchy in this case. In this figure priors have been added. 
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Figure 6: Marginalized octant sensitivity from a combination of NOi^A and T2K, for the case of 
normal and inverted mass hierarchy. The test hierarchy has been left free in this case. In this figure 
priors have been added. 



Figures ^ and ^ depict the octant sensitivity from a combination of NOi^A and T2K 
marginalized over the test parameters and with added priors. The is defined as 

(xL)"""' = min(xNOvA + XT2K + Xprior) (3.6) 

Here (and elsewhere), 'min' denotes a marginalization over the test parameters as outlined 
above. In Figure |5[ the neutrino mass hierarchy is assumed to be known in each case, and 
in Figure ^ the mass hierarchy is taken to be unknown and therefore marginalized over. 
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Table 2: Marginalized octant sensitivity from NOj^A, T2K and a combination of the two experi- 
ments for (5pp — and sin^ 20^3 = 0.1 without added priors (first three columns) and with added 
priors (fourth column). Here normal hierarchy (inverted hierarchy) is assumed. 



These plots are done for four specific true values of the CP phase, 5cp = 0, 7r/2,7r and 
-7r/2. 

Table ^ lists the values of the octant sensitivity from NOz/A and T2K individually 
and from the combination for a specific value 6^p = 0, and a set of true 023 values from 
both octants, without and with added priors, to depict the relative contributions from 
the two experiments. Figure shows how the octant foi^ ^23 ~ varies with the 
test values of 623 and 6c p for the appearance and disappearance channels of NOz^A and 
for T2K. In the figure, ^13 and 6cp are marginalized in the left panel and 613 and 623 
are marginalized in the right panel. Normal hierarchy is assumed in the figure. The 
following features are observed in the sensitivity behaviour of the two experiments and the 
appearance/disappearance channels in each case: 



1. Figure (left panel) shows that the minima for the NOi^A disappearance channel 

^23, 

^23 



occur near ^23'^* = 7r/2 — 6*23, because of the predominant dependence on sin^ 26*2 
while for the appearance channel the minima occur near ^fe*** = 45°, because of the 



sin 023 dependence. In the combination of NOi^A appearance and disappearance 
channels (red solid curve in the figure), the minima are near ^23** = vr/2 — ^231 
following the behaviour of the disappearance channel but the values are enhanced 
due to the contribution from the appearance values at that point. 

The values are asymmetric across the lower and higher octants, as can be seen in 
Table ||. For example, Xmin about 1.5 for NOi/A (appearance + disappearance) 
at (9^^3 = 36° with sin^ 26*^3 = 0.1, Jgp = 0. For 6*^^ = 54°, xLn goes up to about 11 
for NOi^A (app + disapp) for the same value of S^p. 

The values for NOz^A and T2K strongly depends on the true value oi 6cp ■ We 
find for (5^p = tt/2, Xmm for NOi/A (app + disapp) is 5.6 for ^23 = 36° and 7.5 
for ^23 = 54°. This point will be illustrated in more detail later in the context of 



Figures 12 and 13 
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Figure 7: Octant sensitivity from NO^^A appearance and disappearance channels and their com- 
bination and from T2K as a function of test 023 (left panel) and test Sep (right panel). ^13 and 
Sep are marginalized in the left panel and ^13 and 023 are marginalized in the right panel. 



From Figure ^ (right panel) , we see that the ^or the disappearance channel is 
weakly dependent on the test value of Scp, and remains small for all values of test 
5cp. On the other hand, the for appearance channel has a strong dependence 
on test 6cp, which is consistent with the behaviour of the respective probabilities seen 



in Section \2,.2[ In terms of test Sqp, the minima of the combination occur at test Scp 
values close to those for the appearance channel, which is the principal contributor to 
the 6cp dependence of the sensitivity. The increased values of the combined NOi^A 
minima at each test Sqp, can be attributed to the tension between the two channels 
with respect to 9^^. Thus marginalization over dcp reduces the sensitivity to some 
extent because of the effect of 6cp uncertainty in the appearance channel. 

5. The x"^ minima for both appearance and disappearance channels are very low for all 
values of ^23) indicating negligible octant sensitivity from the channels separately. 
This is apparent from the left panel of Figure where the minimum values of for 
individual channels are quite small when marginalized over the entire 623 range. 

6. There is no major tension between NOi^A appearance + disappearance and T2K 
appearance + disappearance. The x^ minima occur at the same values of test sin^ 2^13 
and test 6*23, and for somewhat displaced values of test 6c p- This is because both T2K 
and NOz/A roughly follow the octant behaviour dictated by the vacuum probabilities 
with subleading contributions of the (5cp-dependent terms. The slightly increased 
matter effect contribution for NOz^A does not seem to have a significant effect. The 
combined N0z^A-|-T2K x^ values are nearly equal to the sum of the corresponding 
individual NOi^A and T2K values, as seen in Table For example, Xmin for T2K 
appearance + disappearance is about 3 for 9^^^ = 54°. The corresponding NOi^A 
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value is about 11. For N0z^A+T2K, Xmm about 15 for the same ^23- 

7. The addition of priors to N0i^A+T2K enhances the sensitivity drastically, making 
it (for example) (Xmin)^""'^ ~ ^23 — This underscores the importance of 
precision measurement of ^13 on octant sensitivity. With more precise measurements 
of sin^ 2^13 expected from the reactor experiments, the octant sensitivity is expected 
to get better. In the limit of infinite precision, i.e. if ^13 is held fixed, the goes up 
from 3.9 to 7.7 for (9^^, = 39° and Jcp*'^ = 0. 

8. We see from the Figures ^ and ^ that the is asymmetric between the lower and 
higher octants. The nature of this asymmetry depends on the true 6c p value. For 
instance for 6^p = 0, the values on the LO side are much lower than on the 
HO side, but this is reversed for other 6qp values. Note that the addition of priors 
makes the asymmetry less pronounced since the prior contributions for the lower and 
higher octants are roughly symmetrical. This can be seen by comparing the with and 
without prior values of x^ Table 

3.2 Atmospheric neutrinos 

For our study of octant sensitivity in atmospheric neutrino experiments, we look at the 
following set-ups: 

1. A large magnetized iron detector with an exposure of 500 kT yr, capable of detect- 
ing muon events with charge identification, using the following neutrino energy and 
angular resolution: cte^ = 0.1^/E^, ag^ = 10°. 

Note that such a detector will be constructed by the India-based Neutrino Obser- 
vatory (INO) collaboration [25]. The energy and angular resolutions of muons are 
available from INO simulation code [68,69]. But the work on reconstruction of neu- 
trino energy and angle requiring the resolutions for both muons and hadrons are in 
progress. In our neutrino analysis therefore we use fixed resolutions in terms of neu- 
trino energy and angle as described above. Determination of the octant sensitivity 
of INO using the resolutions obtained from INO simulations is in progress. 

2. A LArTPC with an exposure of 500 kT yr (unless otherwise stated) capable of de- 
tecting muon and electron events. No charge identification is assumed here. The 
angular resolutions are taken to be [43]: 



^e., = 3.2° (3.7) 

For the neutrino energy resolutions, we use the estimated value aE„ = 0.1-v/^. A 
LArTPC has a very good angular resolution since ionization tracks can be trans- 
ported undistorted over distances of several metres in highly purified LAr, allowing 
for excellent direction reconstruction by recording several projective views of the same 
event using wire planes with different orientations [28] . 
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3.3 Octant sensitivity using atmospheric muon events in a magnetized iron 
detector 
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Figure 8: Marginalized octant sensitivity from muon events with a magnetized iron detector (500 
kT yr exposure) , for the case of normal and inverted mass hierarchy. In this figure priors have been 
added. 



The is defined as, 



ixLr'"' = min(xitn, + Xprior) , (3.8) 



where, 

where xitm = xl + xl 

For the atmospheric analysis, the flux and detector systematic uncertainties are included 
using the method of pulls as outlined in [52,70-72]. 

It is well known that the possibility of magnetization gives the iron detectors an excel- 
lent sensitivity to hierarchy [52,68,73-76]. Therefore it is plausible to assume that hierarchy 
would be determined before octant in these detectors. Thus in our analysis of atmospheric 
neutrinos, we assume the neutrino mass hierarchy fixed to be either normal or inverted, 
and do not marginalize over the hierarchy. This is tantamount to the assumption that the 
wrong-octant-wrong hierarchy solutions are excluded. In Figure ^, we present the octant 
sensitivity of a magnetized iron calorimeter detector with an exposure of 500 kT yr. The 
second column of Table |^ displays the values for different values of 623 in the LO. The 
in this case is symmetric about '7r/4, and therefore the values in the higher octant are 
similar. Assuming an inverted hierarchy gives worse results because in this case, matter 
resonance and hence octant sensitivity occurs in the antineutrino component of the event 
spectrum, for which the flux and detection cross-sections are lower. A 2a sensitivity is 
obtained at 623 = 39° if NH is the true hierarchy. 
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3.4 Octant sensitivity using atmospheric events in a LArTPC 

A LArTPC is sensitive to both muon and electron type events and one can study the 
interplay of both type of events in octant sensitivity using atmospheric neutrinos. We give 
the results separately for the atmospheric muon and electron signals as well as for the 
combined analysis. 

In the following analysis we assume that the hierarchy will be determined by the 
combination of reactor, long baseline and INO experiments [68] before a LArTPC can give 
an octant measurement. Therefore we do not marginalize over hierarchy in the wrong 
octant. 

The muon and electron event spectra are related to the oscillation probabilities as 
follows: 

~ {fpfiPfifi + 4>ePefj) X (CC cross — section, exposure, efficiency) 
-^e ~ {4>tJ.Pfie + 4>ePee) X (CC cross — section, exposure, efficiency) (3.9) 

The antineutrino event rates are given by the same expressions with the flux and 
probabilities replaced by their antineutrino counterparts. 
The in this case is given by 

{xlT'°' = min(xita + X^rior) , (3.10) 

where, 

where xLm = xl+fi + xl+e 

The atmospheric muon flux is approximately twice that of the electron flux. Hence 
the behaviour of the muon events is dictated by the muon survival probability and to a 
lesser extent by the Pe/i oscillation probability. The sensitivity derivable from muon events 
is somewhat compromised by the fact that Pe^^ and P^^ shift in opposite directions for a 
transition from one octant to another, as observed in Figure ^. On the other hand, for 
electron events the only source of the octant sensitivity is from P^gj since the electron 
survival probability is independent of 023- Therefore both muon and electron events indi- 
vidually give comparable values of octant sensitivity. The sensitivity from muon events is 
more strongly dependent on the value of ^23) s-nd is therefore higher than the sensitivity 
from electron events for values further away from maximal, getting closer to the latter as 
^23 approaches 45°. This can be seen from Table |3| where we list the values of octant sen- 
sitivity from atmospheric muon and electron events in a LArTPC with an exposure of 500 
kT yr, marginalized over the allowed ranges of the oscillation parameters. The neutrino 
mass hierarchy is fixed to be either normal or inverted. The results presented are for true 
values of ^23 lying in the lower octant. The behaviour for in the upper octant is similar 
and symmetrical about ^23 = ^5°. 

Figure |9| shows the behaviour of the octant sensitivity with the atmospheric neutrino 
signal in a LArTPC for a specific true value (^23)*'^ = 39°, as a function of the test values 
of ^23 7 for both muon (left panel) and electron (right panel) events. For the latter the only 
octant-sensitive contribution is from P^ej which does not suffer from the intrinsic octant 
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Table 3: Marginalized octant sensitivity from atmospheric muon, electron and muon + electron 
events with a LArTPC (500 kT yr exposure), for the case of normal (inverted) mass hierarchy, with 
sin^ 20^3 = 0.1 and 5cp^^ = 0. Priors have not been included here. 

degeneracy at the leading order and the octant sensitivity is due to the sin^ 623 x sin^ 20j^ 
term. Since sin^ 20j^ becomes close to 1 near matter resonance, the octant sensitivity 
increases proportionally with the test value of sin^023- Hence the minimum sensitivity is 
seen in the region near 9^^^^ ~ 45° in the wrong octant. The sohd red curve in the left panel 
is for sin^ 2^13 = for which the contribution comes from P^^ only. The curve exhibits 
no octant sensitivity due to the intrinsic degeneracy coming from the sin^ 2^23 term in 
P^^. This degeneracy is broken for higher values of sin^ 2^13 where both P^^ and Pc^ can 
contribute. As a combined effect of all these factors the octant sensitivity does not have 
a clean proportionality with ^^F*' but can have a minimum anywhere within the range 
45" < 6*2!'^* < (90° — 6*23) depending on the relative weightage of the two terms for specific 
energies, baselines and values of the neutrino parameters. This is reflected by the blue 
dot-dashed curve in the left panel. 

We now study whether combining muon and electron events will give an enhanced 
octant sensitivity compared to that from the muon and electron event spectra individually. 
The last column of Table ^ lists the x^+e values for NH (and IH) . These values are obtained 
by marginalizing over all test parameters, but without adding priors. In Figure |l^ the 
sensitivity is presented with the inclusion of priors. The following points may be noted 
from the listed results: 

1. The values of the octant sensitivity that can be derived from the atmospheric muon 
and electron events separately are relatively low. A combined analysis of muon and 
electron events gives improved sensitivities which are better than the sum of the 
marginalized sensitivities from the two kinds of signals. This is because the muon 
and electron event spectra behave differently as functions of ^23 and hence their 
minima occur at different parameter values. Utilizing this fact, a marginalization 
over the sum of values gives an improved result compared to the individual con- 
tributions. The enhancement in the sensitivity due to this synergy can be anywhere 
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Figure 9: Fixed-parameter octant sensitivity from atmospheric neutrinos with a LArTPC (500 kT 
yr exposure) as a function of test 023 for (sin^ 2^13)4^ = 0.1 and (023)tr — 39°. The left panel is for 
muon events and the right panel is for electron events. 



between 20% and 50%, depending on the true value of 623, as can be seen in Table |. 

The addition of prior information further improves the sensitivity by about 20 - 
40%. The most significant contribution comes from the reduced error range of the 
parameter ^13 from recent reactor results. For an assumed normal hierarchy, a 3a 
signal of the octant may be achieved for = 39°(51°) for a true lower octant (higher 
octant) with this exposure, as seen in Figure ^ 

Assuming an inverted hierarchy gives worse results because in this case, matter reso- 
nance and hence octant sensitivity occurs in the antineutrino component of the event 
spectrum, for which the detection cross-sections are lower. 



In Figure 11, the marginalized octant sensitivity from a LArTPC is plotted as a func- 
tion of detector exposure for two true values of 023- The figure shows that it is possible to 
achieve a 2a sensitivity for = 39° and a 3a sensitivity for ^23 = 36° from a LArTPC 
alone for an exposure as low as 120 kT yr, in the case of a normal hierarchy. 



3.5 Effect of magnetization 



Figure 10 shows the values of octant sensitivity derivable from a LArTPC with and without 
magnetization. The magnetization offers the possibility of charge identification. In the case 
of a LArTPC with magnetization, charge identification details are incorporated as given in 
Refs. [29,30,43], assuming a 100% charge identification capability for muon events and a 
20% charge identification capability in the energy range 1-5 GeV (none for higher energies) 
for electron events. 
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Figure 10: Marginalized octant sensitivity from both muon and electron events with a LArTPC 
(500 kT yr exposure), for the case of normal and inverted mass hierarchy. In this figure priors have 
been added. 
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Figure 11: Octant sensitivity for a non-magnetized LArTPC as a function of exposure, has 
been plotted for normal hierarchy for two values of 023 (true) - 36° and 39° . 



The octant sensitivity without charge identification is seen to be about 80% (60%) of 
that with charge identification capability with a similar exposure, in the case of a normal 
(inverted) mass hierarchy. This difference in behaviour is explained as follows. 

If the true hierarchy is normal, neutrino events have octant sensitivity due to resonant 
matter effects, while antineutrino events do not. Therefore the antineutrino events for both 
octants are almost the same, say k. In the absence of charge-identification, we add these 
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events, and 

2 _ ([jV-f^o + k]- [N^o + k]f 
On the other hand, if we do have charge-identification, then the chisq is simply 

^ NHO ^ k 

In both cases, the numerator is same, but the denominator is more in the former case. 
Therefore charge-identification gives us higher sensitivity. 

Note that when the hierarchy is normal, k comes from the antineutrino events. Because 
of the flux and cross-section being small, is a relatively small number. But if the hierarchy 
is inverted, then k will be because of neutrino events, whose higher flux and cross-section 
will make the denominator quite large. Therefore, the reduction in sensitivity due to loss of 
charge-identification will be more apparent for the inverted hierarchy. Charge identification 
capability has been shown to play a very crucial role in determination of mass hierarchy in 
such detectors [43]. However this does not seem to play such an important role for octant 
determination if we already assume a prior knowledge of the hierarchy. If we assume 
that the hierarchy is not known then the octant sensitivity without charge identification is 
unaffected for 1^23 ~ 45° | < 7°, but for smaller (larger) values of ^23 iii the lower (higher) 
octant, there is a drop of 20-30% in the sensitivity. 

When we compare the octant sensitivity from a magnetized iron detector with the 
LArTPC results without charge identification, the iron detector sensitivities are about 
40% of those with a LArTPC if the hierarchy is normal and about 50% of the LArTPC 
sensitivities if the hierarchy is inverted. The sensitivities from a LArTPC are higher than 
those from a magnetized iron detector . This is because of the possibility of realizing very 
high angular resolutions in the former type of detector and the significant contributions 
from muon as well as electron events. 



3.6 Effect of 5cv 

The dependence of the octant sensitivity on the CP phase 6cp is different for atmospheric 
neutrino experiments and the long baseline experiments considered. Because of the strong 
earth matter effects over a large range of atmospheric neutrino baselines, the behaviour of 
the corresponding oscillation probabilities is governed by the enhanced resonant features 
and their dependence on 5cp is suppressed. On the other hand, for N0i^A/T2K baselines, 
since the matter effects are smaller, 5cp plays a greater role in the probability behaviour 
and hence in the octant sensitivity, because of the degeneracy of the octant with ^13 and 
(5cp as explained previously. 



These features are refiected in Figure 12, in which the marginalized octant sensitivity 



for a LArTPC with atmospheric neutrinos, NOi^A, T2K and their combinations are plotted 
as a function of the true (5cp for ^23 = 39° (LO) and 51° (HO) for both normal and inverted 
mass hierarchies. The 5cp behaviour of the NOz^A and T2K sensitivities are seen to be 
similar, and they follow an opposite behaviour with respect to <5cp for ^23 lying in the lower 
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Figure 12: Marginalized octant sensitivity with priors as a function of true Sqp from NOz^A, T2K, 
NOi^A + T2K and for an atmospheric neutrino experiment with a LArTPC for the case of normal 
and inverted mass hierarchy, for 2 values of in the lower and higher octants. 



and higher octants, i.e. for 6^^ = 39° they are higher in the range 0" < Sqp < 180° and 
for 6*23 = 51° they are higher in the range -180° < Scp^'' < 0° [40]. 

For e^23 = 39° (LO) and NH (the top left panel) the contribution from a LArTPC is 
higher than that from combined N0i^A+T2K excepting for a narrow range 6cp^^ = 90° 
to 135°. For = 51° (HO) and NH (bottom left panel) the same trend is observed, 
except that in this case the combined N0i^A+T2K has a higher sensitivity for the range 
5cp*^ = -15° to -90°. 

For the top (bottom) right panels which are for ^23 = 39° (51°) and IH, the NOz^A -f 
T2K as well as the standalone NOi^A contribution to the sensitivity are higher than the 
atmospheric contribution over almost all values of true 5cp- This is because the sensitivity 



-24- 



Q- 

o 
to 



H On 




1 oU 








9U 




bU 




on 
oU 









-30 




-60 




-90 




-120 




-150 




-180 





37 39 41 43 45 47 49 51 
023 (true) 



o 

CO 



180 
150 
120 
90 
60 
30 

-30 
-60 
-90 
-120 
-150 
-180 



37 39 41 43 45 47 49 51 
(true) 



53 



Figure 13: Wrong octant exclusion for N0j'A+T2K in tlie true 623 - true Sq^ plane. The 
left(right) panel is for NH(IH) as the true hierarchy. The brown(dense-dotted)/red(sparse- 
dotted)/gold(dashed)/green(solid) contours denote 1/2/3/4ct exclusion, respectively. 



from atmospheric neutrinos is less than half for an inverted hierarchy compared to normal 
hierarchy, while for NOz^A and T2K the values are similar. 

The T2K sensitivity is lower than NOz^A in all cases, and lower than the atmospheric 
contribution for most values of true 5cp, besides certain small ranges in the inverted hi- 
erarchy case where it is higher. The combined sensitivity from atmospheric neutrinos and 
NOi^A + T2K follows the true 5cp dependence of the NOz^A + T2K sensitivity. 



Finally, in Figure 13, we show the wrong octant exclusion sensitivity for N0z^A-|-T2K 
in the true ^23 - true 5cp plane. Contours are shown in this plane for 1, 2, 3 and 4(T octant 
exclusion. In the part of the parameter space that is enclosed within a certain contour, it is 
possible to exclude the wrong octant with the corresponding confidence level. As expected, 
the ability of the experiments to exclude the wrong octant is small around true ^23 — 45°. 
As ^23 becomes more non-maximal, the exclusion ability becomes better. The effect of (^cp 
on octant sensitivity for any given true value of ^23 can also be read from these plots. Since 
the results using atmospheric neutrinos are almost independent of JcP) the behaviour of 
the is similar to the plots in Figures |8| and 10 drawn for 5cp = 0- 



3.7 Octant sensitivity from combined analysis of atmospheric electron and 
muon events v^ith NOz^A and T2K 

In this section we present the octant sensitivity from a combined analysis of simulated long 
baseline and atmospheric data. We add x^prior to the combination to take into account 



the future precision measurements on ^13 and 9^^ and minimize the (x 



tot)''™' given as. 



■ / 2 I 2 1 2 I 2 ^ 

mml^XNOvA + Xt2K + XAtm + XpriorJ 



(3.11) 
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Figure 14: Marginalized octant sensitivity from a combination of the atmospheric muon neutrino 
signal in a magnetized iron calorimeter detector (500 kT yr) + NOi^A + T2K, for the case of normal 
and inverted mass hierarchy. In this figure priors have been added. 
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Figure 15: Marginalized octant sensitivity from a combination of the atmospheric muon + electron 
neutrino signal in a LArTPC (500 kT yr, without charge ID) + NO;/A + T2K, for the case of normal 
and inverted mass hierarchy. In this figure priors have been added. 



Figure 14 shows the combined octant sensitivity using the atmospheric muon neutrino sig- 
nal in a magnetized iron detector with NOi^A and T2K. In Figure |l^ the octant sensitivity 
computed from a combination of the atmospheric electron and muon neutrino signal in a 
LArTPC without charge identification capability with NOz^A and T2K is plotted. Finally, 
in Table |^ we list the marginalized octant sensitivity with priors for a LArTPC (500 kT yr 
exposure), a magnetized iron detector (500 kT yr exposure), and a combination of each of 
them with NOz^A + T2K, for the case of both NH and IH. All of the above results are for 
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Table 4: Marginalized octant sensitivity with a LArTPC (500 kT yr exposure), a magnetized iron 
detector (500 kT yr exposure), and a combination of either of them with NOi^A + T2K, for the 
case of normal (inverted) mass hierarchy, with sin^ 20^3 = 0.1 and Scp*^ — 0. Priors have been 
included in this table. 



(5cP*'' = 0. For a magnetized iron detector, the combination with NOi^A + T2K gives a 3a 
sensitivity at 6*23 = 39° for NH. A non-magnetized LArTPC + NOz^A + T2K gives a 4(T 
sensitivity in the same case. There is a tension between the behaviour of the N0z^A/T2K 
octant sensitivity and the octant sensitivity from an atmospheric neutrino experiment as 
a function of test 823, which can be seen by comparing Figures and |. For N0i/A/T2K, 
the minima occur at or close to 9^^^ = tt/2 — 9^^, while for an atmospheric experiment 
the muon events may have minima anywhere between 6^^^ = 45° and it/2 — 6^23 the 
electron events have minima close to 9^^^ = 45°. This synergy leads to an enhancement of 
the octant sensitivity when the NOi^A + T2K and atmospheric neutrino experiments are 
combined. 

The combined values depend strongly on the true value of 5c p following the be- 
haviour of the N0i^A/T2K x^, as can be seen from the green (dot-dashed) curve in Fig- 
ure 12. The LArTPC+N0z^A+T2K combination is seen to reach a sensitivity of > Aa for 
all values of 5cp for NH for the sample values of 6*23 of 39° and 51° considered (the top 
and bottom left panels). For IH and 39° (51°) one achieves close to 2>a (4(t) sensitivity. 
The combination of the magnetized iron detector with N0z^A+T2K is also expected to 
have a similar behaviour with 5c p albeit with a lower sensitivity. The NOi^A and T2K 
contribution play a greater role in this case than that for a LArTPC. The sensitivity for 
N0i^A+T2K is higher for IH than for NH for (Jcp**^ = and ^23 lying in the LO (as seen in 
Table Therefore for these parameter values, the octant sensitivity for a magnetized iron 
detector+NOz^A+T2K combination has comparable values for NH and IH, even though the 
sensitivity from a magnetized iron detector alone is less for IH than for NH. In the case of 
a LArTPC+N0i^A+T2K combination, because of a greater contribution from LArTPC, 
the sensitivity is less for IH than for NH, reflecting the behaviour of the sensitivity from 
LArTPC. This feature can be observed in Table B. 
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4. Summciry and Conclusion 



In this paper we have studied the possibihty of determining the octant of the atmospheric 
mixing angle 623 in the long baseline experiments T2K and NOz^A as well as by atmospheric 
neutrino experiments. While the octant degeneracy conventionally refers to the indistin- 
guishability between 623 and tt/2 ~ 623, this can be generalized to include the whole range 
of allowed value of 623 in the wrong octant, and we consider this generalized definition in 
our analysis. 

We present a probability level discussion on the effect of uncertainty in ^13, Scp and 
values of 623 in the wrong octant for baselines relevant to long baseline and atmospheric 
experiments. Below we summarize the salient features that emerge from our study at the 
probability level: 

• For baselines where matter effects are small, the appearance channel probability 
Pjue displays a degeneracy of the $23 octant with the values of ^13 and 5cp due to 
its dependence on the combination sin^ 623 sin^ 2^13 at leading order as well as its 
subleading c^cp dependence. The disappearance channel probability P^^, suffers from 
an intrinsic octant degeneracy between 623 and 71/2 — 623 due to being a function of 
sin^ 2^23 at leading order. 

• For N0z/A/T2K baselines, there is a strong effect of the uncertainty in 6cf in the 
appearance channel, and of the O23 uncertainty in the disappearance channel. For 
instance, for the NOz/A baseline the degeneracy between 623 = 39° and 51° in P^e is 
lifted only for sin^ 2^13 ^ 0.12 if a variation over the entire range of ^cp is taken into 
account. 

• We find that after including the improved precision in ^13 from the recent reactor 
results, the octant degeneracy with respect to ^13 is largely reduced. Some amount 
of degeneracy with respect to 613 still remains for higher (lower) values of 613 in the 
allowed range for a true lower (higher) octant. Also, in general, information from the 
energy spectrum helps in reducing the degeneracy since for two different energies the 
degeneracy exists for different sets of parameter values. 

• For a sample baseline of 5000 km traversed by the atmospheric neutrinos strong 
resonant matter effects help in lifting both the intrinsic octant degeneracy in P^^ and 
the octant degeneracy with ^13 and Scp in P^e- This is because the term sin^ 29^ 
becomes close to 1 at or near matter resonance, and this makes the leading-order 
term proportional to sin^ ^23 (sin'^^?23) in P/ie (P/i/^) predominate and give distinct 
values in the two octants irrespective of the vacuum value of 613. 

We perform a analysis of the octant sensitivity for T2K/N0z^A and atmospheric 
neutrinos as well as a combined study. For atmospheric neutrinos we consider two types 

of detectors - magnetized iron calorimeter detectors capable of detecting muons and iden- 
tifying their charge and a LArTPC (non-magnetized) which can detect both muons and 
electrons with superior energy and angular resolutions. For N0z/A/T2K we present the 
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results for both cases of known and unknown hierarchy while the atmospheric results are 
assuming hierarchy to be known. The main results are summarized below: 

• The minima for the N0^'A/T2K disappearance channel occur near = ir/2 — 
^23) because of the predominant dependence on sin^ 2^23 5 while for the appearance 
channel the minima occur near ^23'^* = 45°, because of the sin^ ^23 dependence. This 
leads to a tension between the behaviour of the two spectra as a function of ^?2F*, 
so that the combination of appearance and disappearance channels enhances octant 
sensitivity, but even then the is not very high. Combining the NOi^A and T2K 
data leads to a higher octant sensitivity due to the addition of sensitivities from the 
two experiments. What plays a major role in enhancing the octant sensitivity in 
these experiments is the addition of priors, especially on ^13, which helps in ruling 
out the degenerate solutions in the wrong octant. After adding priors one can achieve 
a ~ 2(7 sensitivity at 6*23 = 39° for sin^ 2^13 = 0.1 and 5cp = for both normal and 
inverted hierarchies. 

• A magnetized iron calorimeter gives a 2a sensitivity to the octant for ^?23 = 39° and 
sin^ 2^13 = 0.1 for a true normal hierarchy. A non-magnetized LArTPC can give a 3a 
sensitivity for the same parameter values. The enhanced sensitivity of the LArTPC is 
due to the contribution of the electron events as well as superior energy and angular 
resolutions compared to an iron calorimeter detector. In a LArTPC, combining the 
muon and electron events gives improved sensitivities due to the different behaviour 
of the muon and electron event spectra with respect to 023- If a LArTPC can have 
charge identification capability there will be a 20-40% increase in sensitivity, the 
enhancement being more in the case of an inverted hierarchy than for a normal 
hierarchy. Since we assume hierarchy to be already known, the magnetization of 
LArTPC (which is technically challenging) does not play a significant role. 

• Combining long baseline and atmospheric results can give a significant enhancement 
in the octant sensitivity because of the tension in the behaviour of the N0z^A/T2K 

and atmospheric functions of test ^23) arising due to the different energy and 

baseline ranges involved and the fact that strong matter effects dictate the proba- 
bility behaviour for atmospheric baselines. Since the minima occur at different 
parameter values, there is a synergy between the experiments which gives an in- 
creased octant sensitivity. A combination of NOi/A, T2K and LArTPC (500 kT yr) 
can give a nearly 3a sensitivity for ^23 = 41°, a 4a sensitivity for ^23 = 39° and a 
5a sensitivity for ^23 = 37° for sin^ 2^13 = 0.1 and Scp = in the case of a nor- 
mal hierarchy. The corresponding sensitivities for inverted hierarchy are somewhat 
less. A magnetized iron calorimeter combined with NOz^A and T2K can give ~ 3a 
sensitivity for ^23 = 39° for both NH and IH. 

• The octant sensitivity from an atmospheric neutrino experiment is nearly indepen- 
dent of the true value of 6cp, while the NOi^A and T2K sensitivities arc seen to 
be strongly (5cp-dependent and follow a definite shape as a function of Sqp- The 
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shape gets flipped for a true lower or higher octant. For a normal mass hierarchy, 
the contribution of NOz^A + T2K to the octant sensitivity is lower than that from a 
LArTPC measurement for almost all values of 5cp, while for an inverted hierarchy it 
is higher, due to the worse sensitivity from an atmospheric experiment in the inverted 
hierarchy case. 

In conclusion, we find that the improved precision of ^13 and the different dependence 
on 6*23 in the disappearance and appearance channels and in vacuum and matter probabil- 
ities leads to an enhanced octant sensitivity when long baseline and atmospheric neutrino 
experiments are combined. 
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